Acta Phys. Hung. A 19/1 (2004) 000-000 

HEAVY ION 
PHYSICS 



Heavy Quarks and Heavy Quarkonia as Tests of 
Thermalization 

J.L.Nagle 1 for the PHENIX Collaboration 

1 University of Colorado, 
Boulder, CO 80305, USA 



Received 15 September 2005 

Abstract. We present here a brief summary of the presentation given at the 
"Quark Gluon Plasma Thermalization" Workshop in Vienna, Austria in August 
2005, directly following the International Quark Matter Conference in Hungary. 

Keywords: Quark Gluon Plasma, Heavy Ions, Heavy Flavor 
PACS: 25.75.-q,25.75.Nq 



1. Introduction 

In the PHENIX White Paper [^, we reported the following conclusions: (1) At 
RHIC we have created bulk matter at energy densities well above that predicted by 
lattice QCD for the transition to a Quark Gluon Plasma (QGP). (2) The energy 
density is dominantly equilibrated at very early times (< 2 fm/c), which is when 
the energy density is highest. (3) The bulk matter behaves collectively and as such 
has been described as a nearly perfect fluid. We want to push these conclusions 
further utilizing new data from the large statistics Au + Au and Cu + Cu running 
at RHIC, in particular on heavy quark dynamics and heavy quarkonia suppression 
or lack thereof. 

Almost three years ago, some of us suggested that the PHENIX data on non- 
photonic electrons (presumably from heavy flavor meson decay) may be consistent 
with charm thermalization and hydrodynamic flow At the time, many dismissed 
this hypothesis, and yet now this is the commonly held belief in the field and 
supported by new experimental data. The large charm quark mass means that 
only very strong interactions with high frequency can bring them into equilibrium 
with the light quarks and gluons in the medium. In a calculation by Teaney and 
Moore [[3], they calculate the expected transverse momentum modifications (Raa) 
and momentum anisotropy (1*2) for different charm quark diffusion coefficients - 
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put in as a free parameter in their calculation. Their calculations show that the 
suppression of high transverse momentum charm goes hand in hand with an increase 
in the momentum anisotropy. 

Lattice QCD results show that the confining potential between heavy quarks 
is modified - screened - at high temperatures. At sufficiently high temperatures, 
this screening should suppress bound state formation, such as the J/ip. However, 
recent lattice results indicate that the J/tp spectral function show only modest 
modification near the critical temperature, and thus may not be suppressed until 
significantly higher temperatures. 

2. Experimental Results 

The PHENIX Experiment was designed to measure electrons, muons, photons and 
hadrons utilizing rare event triggers and high data acquisition throughput [0]. For 
the heavy quark and quarkonia results detailed in this proceedings, we utilize the 
particle identification of electrons in two central spectrometers. The acceptance is 
around mid-rapidity —0.35 < rj < +0.35 and electron-pion separation is achieved 
using a Ring Imaging Cerenkov Counter and track matching to an Electromagnetic 
Calorimeter. Also crucial for our measurement is the very low radiation length (< 
0.4%), which keeps the photon conversion background low. In addition, we identify 
muons at forward rapidities 1.2 < \y\ < 2.2 in the PHENIX muon spectrometers 
through a series of interleaved absorbers and active detectors. 

2.1. Open Charm Results 

PHENIX has published results on open charm indicating that the total charm 
yield scales with the number of binary collisions [ 5 . This indicates that charm 
production may be a "hard process" and not suffer large modification due to 
coherence effects. Note that this result does not comment on modification of 
the distribution of charmed hadrons, but only on the scaling of the integrated 
dN/dy(0.5 < pr < 4.0 GeV) near mid-rapidity. In addition, PHENIX has pub- 
lished the first results at RHIC from a modest Au + Au data sample from Run-2 
revealing a non-zero momentum anisotropy (1)2) for non-photonic electrons 

From Run-4, we have collected a significantly larger data sample and have re- 
duced the converter material near the beam-pipe, thus reducing the radiation length 
before our tracking detectors from 1.3% to 0.4%. We use two different methods for 
extracting the non-photonic electron distribution from the initially measured inclu- 
sive electron sample. One method is to subtract off all known photonic contribu- 
tions, using our own measurements of the tt° and r\ as input. The second method is 
making use of a special "converter run" where we purposely increase the radiation 
length around the beam pipe. We can then compare the inclusive electron yield 
between this "converter run" and normal running to determine the photonic con- 
tribution and subtract it away. The first method has larger systematics at low pt 
since the ratio of non-photonic to photonic electrons is small. The second method 



Heavy Quarks and Heavy Quarkonia 



3 



works quite well at low px, but is statistics limited at high px due to the short time 
duration of the "converter run" . Thus, the two methods are quite complementary 
and also agree quite well at intermediate px where both maintain good accuracy. In 
Figurc^Jwe show the PHENIX Preliminary results for non-photonic single electrons 
as a function of px for various Au+Au centrality selections. In addition, the curve is 
the best fit to our non-photonic electron result from proton-proton reactions, scaled 
up by the expected number of binary collisions. Shown in Figure [21 is the nuclear 
modification factor (Raa) for central Au + Au reactions. We observe a significant 
suppression that appears to increase as a function of px and is strongest for the 
most central reactions. 



Au+Au = 200 GeV, Non-Photonic (e + e*)/2 
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Fig. 1. PHENIX Preliminary Au + Au y / s]y]v=200 GeV data for the invariant 
yield of non-photonic electrons versus transverse momentum for various centrality 
selections. The curve is a best fit to the proton-proton yield scaled up by the 
expected number of binary collisions. 



Calculations assuming only radiative charm quark energy loss are able to de- 
scribe the data with varying degrees of success [0 |S] . However, we note that one can 
always arbitrarily increase dN/dy(ghion) or similarly the q value, but then one may 
observe conflicts with light quark/gluon energy loss results or total entropy limits. 
It has been pointed out [ 8 that beauty meson semi-leptonic decay may contribute 
significantly to the non-photonic electrons for px > 3 GeV/c and that the dead- 
cone effect significantly limits bottom quark energy loss. However, we note that for 
bottom quarks, neglecting collisional energy loss, as opposed to radiative, may not 
be well justified. Other calculations include only collisional energy loss [|3J. Better 
experimental constraints on the charm and beauty total cross sections will also be 
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Fig. 2. PHENIX Preliminary central 0-10% Au + Au yfsWN=200 GeV suppression 
factor Raa for non- photonic electrons. 

an important ingredient in understanding where each contributes and at what level. 
PHENIX has a preliminary measurement of the Upsilon in proton-proton reactions 
which is a start at gauging beauty production. 

PHENIX presented preliminary results from Run-4 on the momentum anisotropy 
(V2) for non-photonic electrons, as shown in Figure^ It is now conclusive that these 
electrons have a substantial non-zero anisotropy. We should note that some care is 
warranted in interpreting these results as "charm flow." We believe these electrons 
arc dominated by semi-leptonic decay of charm mesons and an additional contribu- 
tion from the decay of beauty mesons at higher px- Due to the decay kinematics, 
the electron carries only a fraction of the meson px and the effective vi is reduced 
for low px electrons as the A(f> between the parent and daughter particle is effec- 
tively like an additional reaction plane smearing contribution. We note that the 
data gives an indication for a decrease in t> 2 above a px ~ 2 GeV/c. This could 
be the result of a decrease in charm quark flow as predicted in [ |3] or due to the 
emergence of beauty contributions. 



2.2. Quarkonia Results 

In studying closed charm or beauty in heavy ion reactions, we are interested in the 
interaction between the heavy quark and antiquark and the surrounding medium. 
However, we must also keep in mind that although the state likely begins as par- 
tons, it must transform itself into a hadron before being directly observed. The 
promise of insight from quarkonia measurements is tempered by the large number 
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Fig. 3. PHENIX Run-2 published and Run-4 preliminary t>2 for non-photonic 
electrons as measured in minimum bias Au + Au 200 GeV reactions. 



of possible effects impacting their final yield. There may be nuclear modification to 
the incoming parton distributions (e.g. shadowing, anti-shadowing, EMC,...) that 
may impact the exact scaling of total heavy flavor production. After initial creation 
of the quark antiquark pair, they are bombarded by the "back-side" of the two 
nuclei. These nucleons and the quarks and gluons inside them may break up the 
heavy cc pair - a process referred to as normal nuclear absorption. Then the cc pair 
can interact with the surrounding medium either at the partonic or hadronic level. 

The PHENIX experiment has submitted for publication results on J/ip produc- 
tion in proton-proton and deuteron-Aw collisions at RHIC energies over a broad 
range in rapidity —2.2 < y < +2.2 [12]. In comparing our proton-proton and 
deuteron-Au results, we find a very modest suppression (or none within errors) of 
J/tp at mid-rapidity relative to binary scaling and a possible larger suppression (of 
order 20%) at forward rapidity giving a hint of gluon shadowing effects. Future 
higher statistics deuteron-Au data will be required for more precise conclusions, 
but the current data are consistent with a J/T/>(precursor)-nucleon breakup cross 
section of order 1-3 mb. 

Previously, PHENIX had published only a very low statistics result on J/tp 
production in heavy ion reactions [E3- From the Run-4 and Run-5 high statistics 
samples, the PHENIX experiment has presented preliminary results on J/ip pro- 
duction in Au + Au and Cu + Cu reactions as a function of collision centrality, 
transverse momentum and rapidity. All results are shown together in Figure ^ in 
terms of the nuclear modification factor Raa as a function of the number of partic- 
ipating nucleons. We observe a suppression of J/tjj yields relative to binary scaling. 
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The suppression measured at mid-rapidity (via the dielectron channel) is compa- 
rable within statistical and systematic errors of that measured at forward rapidity 
(via the dimuon channel). 

We overlay the PHENIX data with three different theoretical and experimental 
comparisons. First, shown as the upper two red curves are calculations assuming 
only nuclear modification of parton distribution functions (EKS98) and normal 
nuclear absorption with a cross section a = 3 mb, which is at the limit of agreement 
with our deuteron-Au data [E]. These calculations for y = and y = 2 appear to 
underpredict the level of suppression for the more central Cu + Cu and Au + Au 
data. Next we show various calculations assuming further suppression of the J/ip 
due to comover absorption or disassociation due to screening [ E| ■ These three 
particular calculations were all matched to the J/ip suppression observed at lower 
energies, and since the RHIC energy density is a factor of 2-3 higher, they all 
substantially overpredict the level of suppression. Finally, we show the experimental 
data from lower energy from experiment NA50, normalized to yield Raa = 1 for 
the most peripheral events [Ej- Although one must take seriously the current 
PHENIX systematic errors, the general agreement with the lower energy centrality 
dependence is striking. There are a variety of theoretical calculations invoking J/ip 
regeneration or coalescence [I14|. These calculations qualitatively predict very large 
initial suppression of J/ip, which is compensated for by later re-formation. These 
calculations may give a better description of the data, but must be checked in the 
context of the J/ip proton- proton cross section and the input charm quark cross 
section and distributions [I15|. 

We have also presented PHENIX preliminary transverse momentum and ra- 
pidity distributions of J/ip from Au + Au and Cu + Cu reactions. We find no 
large modifications of these distributions in comparing proton-proton to heavy ion 
reactions, but further quantification of this conclusion requires pushing down our 
current systematic errors. It is notable that since charm quarks show a suppression 
at high pt, one might expect J/ip regeneration to lead to a significant distortion of 
the transverse momentum distribution. Predictions of narrower rapidity distribu- 
tions have also been made. We show in Figure [5] Raa in 17 bins in centrality from 
our Cu + Cu data set, which emphasizes that we have excellent statistics to explore 
various dependencies. For example, it was suggested that perhaps only the \c is 
suppressed and that is why the NA50 and PHENIX suppression patterns appear 
similar. If this were the case one might naively expect the suppression onset to 
occur at RHIC energies in mid-central Cu + Cu, whereas our data show a slow 
onset of the full suppression seen for the most central events. Whether the J/ip is 
really suppressed in medium and then regenerated or is never much suppressed at 
all is a question we hope to answer with the reduction of our current statistical and 
systematic errors and future measurements. 
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Fig. 4. PHENIX Preliminary Au + Au and Cu + Cu 200 GeV nuclear modification 
Raa for J/ip. Various theoretical predictions and experimental data as described 
in the text are shown compared with the PHENIX data. 
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Fig. 5. PHENIX Preliminary R AA for J/tp in Cu + Cu 200 GeV reactions. 
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3. Conclusions 

In summary, there is a wealth of new PHENIX data on heavy quarks and heavy 
quarkonia. We will work hard to push these results to submitted publications. 
Charm is a very optimal probe of thermalization and properties of the medium, but 
the price for this may well be the loss of a probe via quarkonia for deconfinement. 
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